Rod-like Polyelectrolytes in Presence of Monovalent Salt 



Paulo S. Kuhn, Yan Levin* , and Marcia C. Barbosa 
Instituto de Fisica, Universidade Federal do Rio Grande do Sul 
Caixa Postal 15051, CEP 91501-970 , Porto Alegre, RS, Brazil 

We investigate the properties of rigid polyelectrolyte solutions in presence of monovalent salt. 
The free energy within the Debye-Hiickel-Bjerrum (DHBj) theory [M. E. Fisher and Y. Levin, Phys. 
Rev. Lett. 71, 3826 (1993)] is constructed. It is found that at thermodynamic equilibrium the 
polyelectrolyte solution consists of clusters composed of one polyion and various counterions. The 
distribution of the cluster densities is determined by finding the minimum of the Helmholtz free 
energy. The osmotic pressure and the average charge of the cluster are found and their dependence 
on Manning parameter £ is elucidated. A good agreement with the experimental results is obtained. 



1. INTRODUCTION 

Over the last three decades polyelectrolyte solutions 
have found a number of practical applications ranging 
from super-absorbents and hair conditioners to water 
treatment. Many biologically important molecules such 
as DNA, are also polyelectrolytes. Not withstanding 
their practical importance our theoretical understanding 
of the behavior of these complex molecules is still quite 
rudimentary [|l| . The fundamental problem which makes 
the study of polyelectrolytes so much more difficult than 
that of regular polymers is the long-ranged nature of the 
Coulomb force. The scaling theories, which have proven 
so useful for simple polymers have, so far, failed in the 
case of polyelectrolyte solutions. Even the mean-field 
theories are extremely hard to construct, as one tries to 
take a realistic account of the long-ranged electrostatic 
interactions. 

In general a polyion is a polymer some of whose 
monomers are ionized (all of the ionized monomers have 
the same sign of charge) . In most practical applications 
the polyions are dissolved in some solvent, usually wa- 
ter, causing a strong interaction between the charged 
monomers and the counterions. Besides this, already 
very complex interaction, in the case of a flexible chain 
one must also take into account the conformational de- 
grees of freedom of the polymer, as well as the inter- 
actions between the macromolecules inside the solution. 
The full problem is extremely difficult to study, however, 
some simulations have been attempted [||. A somewhat 
simpler problem, an answer to which we shall attempt to 
elucidate in this paper, is the properties of rigid polyelec- 
trolytes. An example of rigid polyelectrolyte is a solution 
of DNA segments. A major simplification resulting from 
restricting our attention to rigid polyions is that these 
can be modeled as cylinders, thus allowing us to bypass 



the complication of taking a full account of the confor- 
mational degrees of freedom. At high densities the rigid 
molecules tend to align forming a smectic phase. The 
periodic structure of the smetic phase allows us a major 
simplification of studying one polyion inside a Wigner- 
Seitz cell, since all of the cells are identical || . At low vol- 
ume fractions, when the solution is isotropic, the Wigner- 
Seitz picture is no longer valid, and a new approach must 
be found. On the other hand, it is well known that at infi- 
nite dilution the polyelectrolyte solutions obey the Man- 
ning limiting laws Q. Although derived in a somewhat 
ad hoc manner, the limiting laws have proven extremely 
successful and provide a "boundary condition" that any 
theory of polyelectrolytes must satisfy. Recently we have 
extended the Debye-Huckel-Bjcrrum (DHBj) theory (|] 
to rigid polyelectrolyte solutions and found that at infi- 
nite dilution it, indeed, reduces to the Manning limiting 
laws j6|. However, consideri ng the higher density cor- 
rections, we have found that the limiting laws should ap- 
ply only at extremely low densities, since the corrections 
scale as l/lnp, where p is the density of polyions [||. 
What can then account for the success of limiting laws 
at densities which are not so small? In this paper we 
shall extend the DHBj theory to realistic concentrations 
encountered in most applications. 

DHBj theory has proven to be successful in a vari- 
ety of systems whose dominant interactions are due to 
the long-ranged Coulomb potential The fundamen- 
tal idea behind the theory is that the nonlinearities which 
are omitted in the process of linearization of the Poisson- 
Boltzmann equation can be re-introduced into the the- 
ory through a thermodynamic assumption that the op- 
positely charged particles can associate, forming clusters. 
The density of these clusters will, then, be determined by 
the condition that at equilibrium the total free energy of 
the system must be minimum. In this sense, the DHBj 
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theory provides a kind of variational approximation to 
the complete, yet unknown, theory. By virtue of being 
linear, DHBj theory bypasses the internal inconsistencies 
which complicate many of the nonlinear theories of ionic 
solutions ||. 

In the case of symmetric electrolytes, the clusters cor- 
respond to the dipolar pairs formed when two oppo- 
sitely charged ions come into a close contact 0. These 
are exactly the kind of configurations which become 
"undervalued" when one linearizes the Boltzmann fac- 
tor. For the polyelectrolyte, the polyions of which can 
carry a charge which can be many thousands units 
(unit=electron charge), there can be a large variety of 
clusters. In general we expect that there will be clusters 
consisting of one polyion and one counterion, one polyion 
and two counterions, one polyion and three counterions, 
etc. The successful theory should allow us to calculate 
the full distribution of cluster densities. In the previous 
work we have demonstrated that in the limit of infinite 
dilution only one characteristic cluster size is thermody- 
namically stable |6|. This is the same as was previously 
postulated by Manning [Q. However, as one moves away 
from infinite dilution, one must consider the full distribu 
tion of cluster sizes ||, and this will be the objective of 
this paper. 



2. THE MODEL 



We shall first define the Primitive Model of Polyelec- 
trolyte (PMP) ||. Our system will consist of long cylin- 
dric polyions, the spherical counterions, and salt, inside 
a volume V. The polyions, of length L and density p p , 
have diameter a p , and charge — Zq uniformly distributed 
along the length of the cylinder. The distance between 
charged groups is b = L/Z. The counterions will be mod- 
eled as rigid spheres of diameter a c and charge +q located 
at the center. The overall charge neutrality requires that 
the density of counterions be Zp p . In addition to the 
polyions and the counterions we shall allow for the pres- 
ence of salt of density p sa u ■ In this work we shall consider 
a simple monovalent salt and, in particular, treat the pos- 
itive ions of salt as identical to the counterions. We shall, 
therefore, consider only two species of small ions, the ones 
with charge +q, which we shall denote counterions indis- 
criminantly of whether they are derived from the p oly- 
mer or the salt, and the coions with charge — q, derived 
from dissociation of salt molecules. The solvent will be 
modeled as a continuum medium of dielectric constant D. 
The distance of closest approach, a, between a polyion 
and a counterion is a = (a p + a c )/2. This represents an 
"exclusion cylinder" due to the hard-core repulsion (see 
Fig.l). 
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FIG. 1. A polyion of (cylindric) diameter a p and length 
i>a surrounded by spherical counterions and coions of di- 
ameter a c . The charge spacing is b = L/Z, and the radius of 
the exclusion cylinder is a = (a p + a c )/2. 

The strong electrostatic attraction between the 
polyions and the counterions will result in some counteri- 
ons becoming associated with the polyion molecules, pro- 
ducing clusters consisting of one polyion molecule and n 
counterions. Just as in the process of micellization, that 
is commonly observed in amphiphilic systems, we expect 
that there will be a distribution of cluster sizes. Thus, in- 
side the polyelectrolyte solution, we shall encounter some 
free polyions with no attached counterions. We denote 
these zero-clusters of density po . We shall also encounter 
clusters consisting of one polyion and one counterion, 
one polyion and two counterions, etc. The density of n- 
clustcrs is p n , with n ranging from 1 to Z. Since not all 
of the counterions condense onto the polyions, for any 
non-zero temperature some free, unassociated, counteri- 
ons remain in the solution. We denote the density of 
unassociated counterions p + . The conservation of the 
total number of par tides leads to two equations, namely 
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and 
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P+ = Psalt + Zp p - 22 n Pn 
n=0 



(2) 



Since the coions do not participate in the association, 
their density remains unchanged, p_ = p sa it- 

The goal of the theory is, therefore, to determine the 
distribution of cluster sizes. The complete thermody- 
namic information about the system is contained in the 
Helmholtz free energy. The condition that the free en- 
ergy must be minimum will allow us to determine the 
distribution of cluster densities. Once the distribution 
is ascertained, all the thermodynamic functions of the 
system can be found through the appropriate operations 
on the free energy. For example, the pressure inside the 
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FIG. 2. Two infinite polyions at distance p and with rela- 
tive orientation 9. The distance between the line elements i 



and j is r = (i — j cos ( 



f+p 2 + (j sine) 2 . 



The eqs. (7) and (8) reduce to a one-dimensional prob- 
lem, since they are symmetric in the angle ip, and the z 
dependence is suppressed due to the large extent of the 
polyion. The only cylindrical variable left is the radial 
distance r from the polyion. The appropriate boundary 
conditions are: the vanishing of the potential for large 
values of r and continuity of the potential and the elec- 
tric field at r = a. With these conditions we find [0] 

€ = -yM^) + yeWifoW, r<a, (9) 



^ t = ^Q(m)K (Kr), r>a, 



(10) 
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FIG. 3. A sample of charge distribution on a polyion. Each 
site contains 1 or counterions associated to it, so the charge 
of the site is or —q, respectively. 



where 
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is the second-class modified Bcsscl function of 
n. The geometric factor 0(«a) arises from the 
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order 

absence of screening inside r < a fll2fl . If instead of a 
cylinder a polyion would be a line of charge, i.e., a = 0, 
then the geometrical factor would become 0(0) = 1. 

The first term in eq. (9) is the potential of an infinite 
line of charge, while the second is the potential of the 
polyion due to the presence of other polyions and coun- 
terions. The electrostatic energy of a particular cluster 
of size n may now be calculated through 
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where A$< n > = <S> { ™ ] + ^ln(r/o), for r < a; A$M = 
$i" ] t , for r > a. That is, we subtract the logarithmic 
potential produced by the line of charge, since it only 
contributes to the self energy of the cluster. The charge 
density is given by eqs. (5) and (6) in linearized form. 
Evaluating the integral, we find 11311 
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The electrostatic free energy is found through the 
Debye charging process, in which all the particles are 
charged simultaneously from to their final charge 
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The free energy density due to the polyion-counterion 
interaction becomes 
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with 



f Ka xK (x) 2 . . 

I(na) = I dx R ^ 2 . (16) 



inclination between the lines, and p is the perpendicular 
distance of separation between the molecules (see Fig. 2).. 
Note that to simplify the calculations we have extended 
the limits of integration to infinity Jlq| . Changing the 
variables, 



4. THE COUNTERION-COION AND THE 
POLYION-POLYION INTERACTIONS 

4.1 The counterion-coion interaction 

In the previous calculation of the polyion-counterion 
free energy we have neglected the correlational effects 
arising due to the counterion-counterion, counterion- 
coion, and the coion-coion interactions. To leading or- 
der the contribution to the excess free energy arising 
from the correlations in the charge density distribution 
can be calculated using the usual Debye-Hiickel theory. 
To this end we consider the counterions and the coions 
moving in a neutralizing background produced by the 
free polyions and clusters. It is then a simple matter to 
demonstrate that the correlational free energy takes on 
the usual Debye-Hiickel form fl(3|f| 
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4.2 The polyion-polyion interaction 

The long-ranged Coulomb interaction between the 
macroions inside the polyelectrolyte solution will be 
screened by the mobile counterions and coions, produc- 
ing a short-ranged effective potential. Following Onsager 
||15f , the energy of interaction between two lines of charge 
in an ionic sea is obtained by integrating the interaction 
of elements of charge, which is given by a Yukawa poten- 
tial, 
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Here o~ m is the linear charge density of a polyion with m 
counterions associated, k 2 = Airpiq 2 j3/D is the inverse 
Debye length, eq. (8), di is the line element of a polyion, 
and r is the distance between the line elements, 



r 2 = (i - j cos 0) 2 + p 2 + (j sin ( 
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Here i and j denote the distance of a line element to the 
center of each molecule, < 9 < it is the relative angle of 
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we have 



: p 2 + S 2 + t 2 
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Hence, in terms of s and t the potential becomes 
c m cr„ f ds dt exp [-n(p 2 + s 2 + t 2 ) 1/2 ] 
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The integral can be performed easily by going to polar 
coordinates, and we find fL5| 
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So far in this section we have treated the molecules 
as lines of charge. A more realistic approximation is to 
model them as cylinders of radius a. As was mentioned 
following eq. (10), one can go from the line of charge to 
the cylinder by replacing o~ n — > a n O(na). In the case of 
spherical particles this transformation leads to the DLVO 
potential of interaction between two macroions Jl2| . We 
shall conjecture that the same procedure remains valid 
for cylinders at asymptotically large separations. With 
this geometrical factor, the pair potential of interaction 
between two clusters becomes 
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It is important to note that this is an asymptotic po- 
tential, valid for p a. When two polyions are in a 
close proximity of each other, the fluctuations in their 
condensed layers can become correlated, producing an 
effective attractive interaction 16 1. Though extremely 



interesting, this short-distance effect should not have any 
significant influence on the thermodynamic properties of 
the polyelectrolyte solution determined mostly by the 
counterions and their interactions with the polyions (see 
Fig.7). 

Since the effective potential is short ranged, we can 
account for the polyion-polyion contribution to the free 
energy using a Van der Waals type of approach. To this 
end the polyion-polyion contribution is expressed as a 
second virial term averaged over the relative angle sus- 
tained by two molecules, 
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The result is 
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Substituting (0) we obtain 
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FIG. 4. Distribution m(s) of associated counterions on a n-cluster for Z = 100, and several values of n. (a) n = 10; (b) 
n = 40; (c) n = 70; (d) n = 90. 



5. THE MIXING FREE ENERGY 

The entropic contribution to the free energy can be 
calculated using the Flory theory Jill . Thus the increase 
in free energy due to mixing of various species is the sum 
of their ideal free energies, 



Pf ent = Y,\pt-ptfoU>t/Ct 



The volume fractions are 
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for the n-clusters, the counterions, and the coions respec- 
tively; Ct is the internal partition function of an isolated 



specie t. For structureless particles, 
while for clusters 
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The Hamiltonian H represents the electrostatic energy 
of interaction between the n condensed counterions lo- 
cated at sites Sj, i — 1, n, and the negatively charged 
monomers of the polyion (see Fig. 3), 
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TABLE I. Free energy of association for Z — 100 and 
£ = 2.283. F R [fh(s)] is the free energy obtained using the 
variational counterion profile. F R [m] is the free energy with 
a uniform profile m = n/Z. 





n 


F R \m(s)} 

L \ / J 


F R \m] 

L J 


10 


— 91 « 44^800 


— 91 4 1 441 79 


20 


-399.827094 


-394.192426 


30 


-554.221516 


-548.635357 


40 


-683.959854 


-679.127270 


50 


-790.048546 


-786.298434 


60 


-872.900058 


-870.322928 


70 


-932.493432 


-931.026672 


80 


-968.353696 


-967.779398 


90 


-979.023957 


-978.926802 





We are labeling the Z sites on the polyion with the 
index s, so 1 < s < Z; q(s) is the net charge on site s, 
which is — q (no counterion associated) or (one counte- 
rion associated); r(s) is the distance (in units of b) of site 
s from one of the ends. We assume that the only effect 
of a condensed counterion is to neutralize the monomer 
charge. Furthermore, the electrostatic repulsion between 
the counterions will prevent that a site is occupied by 
more than one counterion. H[—q] is the energy of the 
reference state in which no counterions are condensed, 
and each site has the charge — q. 

At this point we need to evaluate the internal par- 
tition function. Unfortunately even this comparatively 
simple problem is very difficult to solve due to the long- 
ranged nature of the Coulomb force. In the absence of 
the exact solution we shall resort to a mean-field type of 
approach. We first define the free energy of association 
as —[3F R = ln( n . It is easy to show that this energy is 
bounded from above, 

F R < F R = F R + (H- H ) Q - H[-q] , (34) 

where H is an arbitrary Hamiltonian and F R is the free 
energy corresponding to H - This is a well known Gibbs- 
Bogoliubov bound on free energy. We shall choose Ho to 
be of one-body form so as to facilitate the calculations, 

#o [?(*)] = ^( s M s ) . ?( s ) = . ( 35 ) 

s 

where tp(s) is the effective position dependent mean field 
to be determined. The partition function corresponding 
to Hq is 
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and the free energy, 

F o fl = -^E ln { ex p[^(s)] + i} 

The average charge on site s is 

qexp[/3qip(s)} 



^From this we obtain 
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exp [f3qtp{s)] = - 



P( s ) 
p(s) + q 



Therefore, 
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where m(s) = + 1 is the average number of asso- 
ciated counterions on site s. 

To find the optimal bound on free energy we shall min- 
imize it with respect to m(s) ) i.e., 5F R / 'Sm(s) = 0, un- 
der the constraint J2 S m i s ) = n - The resulting distribu- 
tion m(s), which is obtained by numerically solving the 
minimization problem (see Fig. 4), corresponds (within 
mean field theory) to the equilibrium density profile, i.e., 
(q(s))/q+ 1 w m(s). Without the knowledge of the exact 
solution we shall approximate the free energy of associa- 
tion by its optimal bound, F R ss F [fh(s)]. In practice, 
however, even this simplified procedure will make finding 
the distribution of cluster sizes very slow, since for each 
cluster size n we will need to find m(s) by numerically 
solving the minimization problem. Instead, we noted 
that outside the endcap regions, m(s) is more or less uni- 
form and can be approximated by m(s) « m = n/Z (see 



Fig. 4). With this approximation the sums in eq. (40) 
can be carried out explicitly and we find (£ = (3q 2 /Db is 
the Manning parameter) 

/3F R [m] = £ (m 2 - 2m) {Z [ip(Z) - V(l)] -Z+l} 

+ Z[(l -m.)ln(l -to) + to In. to] , (41) 

where ip(n) is the digamma function fl8|| . We have com- 
pared the free energy obtained using this uniform coun- 
terion profile with that obtained using the exact solution 
of the variational problem and found that they agree to 
within a few percent (see Table 1). We shall, therefore, 
approximate the internal partition function of a rt-cluster 
as ( n — exp(—/3F R [m]) with m = n/Z. Inserting this 
into eq. (28) completes the calculation of the entropic 
contribution to the free energy inside the polyelectrolyte 
solution. 
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FIG. 5. (a) Cluster size distribution p n for polyion density p p = 10 _4 M, salt p sa it = 10 _3 M, and Z — 100; the polyion and 
counterion diameters are 12A and 9.3A, respectively. The average number of associated counterions is n* = 72.173, while the 
Manning infinite-dilution result is um = 56.198. (b) Cluster size distributions p n for two different densities of polyions. One 
should note the strong insensitivity of the distribution to the variation in density of polyions. Note changing scale. 



6. THE HELMHOLTZ FREE ENERGY 

The full Hclmholtz free energy for the polyelectrolyte 
plus salt system is 

j _ jent _|_ jpc _|_ jee _|_ jpp ^2~j 

^From the free energy all the thermodynamic proper- 
ties of the solution can, in principle, be found. Mini- 
mization of the free energy under the constraint of fixed 
number of particles, eqs. ([!]) and (||), reduces to the 
law of mass action, uq + nu + = u n . The chemical po- 
tentials of all species are obtained from the free energy, 



fit = —df/dpt, and in thermodynamic equilibrium we 
find that 

K = ^0 4>l Cn exp[/?(^ + - tf)] , (43) 

where the excess chemical potential is uf x = — df ex /dp t ; 
jex _ jpc _|_ jee _|_ pp ^ an( j |_ ne are ^ e volume frac- 
tions defined earlier. The above expression (43) consists 
of Z + 1 equations for the Z + l densities involved, and 
its solution represents a formidable task. However, the 
numerical treatment of this problem is feasible, and em- 
ploying an iterative algorithm we were able to find the 
distribution of cluster sizes (see Fig. 5). 
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With the cluster densities determined the osmotic pres- 
sure may be calculated, and the results for several values 
of densities of salt are plotted in Fig. 6. A good agreement 
with Manning's low-density limit and the experimental 
results of Alexandrowicz p9[ is obtained. The relative 
contribution of each term of the free energy can be seen 
in Fig. 7. Note that the polyion-polyion contribution to 
the osmotic pressure is negative, what can be interpreted 
as an effective induced attraction between the polyions. 
The characteristic size of a cluster is n* = J2n n Pn/pp- 
In Fig. 8 we demonstrate the dependence of the average 



charge of a cluster, 

on and compare it to the infinite-dilution counterion 
condensation theory B. We note that the sharp coun- 
terion condensation transition disappears at finite den- 
sities, replaced by a smooth crossover. Finally, in Fig. 9 
we demonstrate the relative insensitivity of the average 
cluster charge on the density of added salt. 
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FIG. 6. Osmotic pressure for a polyelectrolyte solution (solid line) for three different concentrations of salt: (a)p aa » = 10 _3 M; 
(b)p sa it = 1CP 2 M; {c)p sa it = 10 — 1 M. We have defined X = Zp p /p sa it- The characteristics of the polyions are: Z = 100, 
£ = 2.283, a p — 12A, a c — 9.3A, a = 10.65A. The dashed line is the infinite-dilution limit of Manning, while the squares, 
triangles and circles are the experimental results of Alexandrowicz (see text for the references). 
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7. CONCLUSIONS 



We have constructed a DHBj theory for an isotropic 
polyelectrolyte solution. We find that, at equilibrium, 
the solution consists of free polyions, free counterions, 
and clusters composed of one polyion and 1 < n < Z 
counterions. The distribution of cluster sizes was cal- 
culated explicitly as was the osmotic pressure. A good 
agreement with the experimental data of Alexandrow- 
icz was obtained. Finally, it is important to note that 
the sharp counterion condensation observed at £ = 1 
for infinite dilution, gets "smoothed out" for finite den- 
sities. This is indeed what one might have expected 
based on the modern theory of phase transitions. It is 
well known that the higher-order phase transitions (sec- 
ond and above) are associated with a diverging length 
scale. In the ionic systems the length scale relevant to the 
counterion condensation phenomena is the Debye screen- 
ing length. In the case of polyelectrolyte solutions, this 
length will stay finite for any non-zero temperature. In 
this respect the counter ion condensation is very different 
from the Kosterlitz-Thouless (KT) metal-insulator tran- 
sition observed in two-dimensional plasma |pcj| . The KT 
transition is a real thermodynamic phase transition at 
which the Debye length becomes infinite. It is important 
to stress this difference in view of the recent specula- 
tions that the counterion condensation is a KT-like phase 
transition Instead we should compare the counte- 

rion association with a micellar formation in amphiphilic 
systems. 
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FIG. 7. Relative contributions to the osmotic pressure 
inside a polyelectrolyte solution (p sa it = 10 _3 M, see 
caption to figure 6 for details) in order of importance: 
mixing (dot-dashed line); polyion-counterion (dotted line); 
polyion-polyion (dashed line); counterion-coion (solid line). 
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FIG. 8. Dependence of the effective charge density a e f / /ao 
on the inverse Manning parameter l/£. Two densities of 
polyions are considered: (a) p p — 10 _4 M; (b) p p — 1CP 14 M. 
The dashed line is the Manning's infinite dilution theory. 
Note that at infinite dilution there is no counterion associ- 
ation for £ < 1, and the transition appears sharply at £ = 1, 
with a e ff /ao = 1 — l/£ (see text in Conclusion). 
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FIG. 9. Dependence of the relative effective charge density 
<r e ff /(To on the density of salt p aa it- The density of polyion 
is fixed: p v = 10" 8 M. 
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